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A method for renormalization of the Casimir energy of confined fermion fields in (1 + \)D 
is proposed. It is based on the extraction of singularities which appear as poles at the point 
of physical value of the regularization parameter, and subsequent compensation of them by 
means of redefinition of the "bare" constants. A finite ground state energy of the two-phase 
hybrid fermion bag model with chiral boson-fermion interaction is calculated as the function 
of the bag's size. 



1 Introduction 

Nowadays a number of methods is developed for evaluation of the Casimir energy for sys- 
tems with quantized fields under nontrivial boundary conditions [fll. The most known of 
them are the Green function method, the zeta function approach, the application of contour 
integration, the multiple scattering method, and the direct mode summation with thermal 
regularization j| (for recent discussion see 0). An important field of application for these 
techniques is an investigation of ground state energy for the quark bag models with account 
of one- loop corrections due to the filled Dirac's negative energy sea Q. In this paper we 
propose an approach for explicit calculation of ground state energy for the models with con- 
fined fermion field coupled to the boson (scalar) field in a certain spatial region, so called 
hybrid bag models ||. This method allows to extract the singular terms from the divergent 
sums of eigenvalues u n and obtain the finite Casimir energy of fermionic sea as the result 
of absorbing singularities into the "dressed" model parameters. The realization of a such 
strategy requires, as a rule, an including of contact terms into the "classical" expression for 
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the energy of system. The form of these terms is determined by their dependence on the 
geometric parameters ||. By means of this method, we renormalize the ground state energy 
of the hybrid chiral bag model in (1 + 1)D and study it as a function of the bag's size. 
It's shown that this function has unique minimum what means that there exists a stable 
configuration for the considered model. 

2 Regular izat ion 

The ground state energy of fermion field is defined as the vacuum expectation value of the 
Hamiltonian: 

(0\H\0) = E = -| ( J2 "n - $>*)■ (!) 

Provided that the spectrum uj n possesses the symmetry u> n — > —u> n , the vacuum energy is 

E = - J2 w«- (2) 

UJn>0 

This sum diverges. On of the ways to regularize it is to use the exponential cutoff J7|, |8|] what 
yields 

E^{r) = lim ( - £ <«*>-™») = lim ± £ e"~», (3) 



where r = f/i™ 1 and /x has the dimension of mass. This regularization can be applied for 
calculation of the Casimir energy in several simple situations |l|, [F|, ||, but usually only for 
systems with an explicitly known spectrum u n . In this paper we'll try to show that one can 
succeed in calculations even if this is not the case. Suppose that an asymptotical expansion 
of uj n for n> 1 | can be found for a model with an unknown spectrum. This expansion 
reads 

— oo ry / 1 \ 

£ £W = fiin + n + — - + — ) . (4) 



— oo 



i=1 n \n< 



These are the three leading terms in (4) which determine the divergencies in (2). To begin, 
consider the case when the expansion (4) contains only two terms (what really takes place, 
e. g. for the free massless fermion field in (l+l)-bag): 

w? = n in + n . (5) 

Then the regularized energy is J7J 

EF{t) = lim A j: e-^ +n °) - wj, = lim ± (e^a^r)) - co , (6) 

n=l e - 1 - 

The term ujq in (6) is written separately because the expansion (4) is not valid for it. By 
virtue of 

_L_ = y E± x k-i (8) 



where B^ are the Bernoulli numbers, one gets 

^« = -^ + l + T + ^-»*- (9) 

This expression contains the quadratic divergence E quad {r) = —72^- = — ^T depending on 
the arbitrary mass \x and the geometric parameters of the bag from f^. 
Consider now the case when 

4 2) = fim + fi + ^ zi • (10) 

n 
In this situation the regularized energy is 

E? p (T) = \imJ-a 2 (T)-co , (11) 

a 2 (r) = J2 e" T ( nin+f7(,+ ^) = e" Tn ° £ g^^e"^- 1 /" . (12) 

n=l n=l 

One easily see that only the two leading terms in the expansion of Q- Tn -^/ n in powers of r 
yield a no n- vanishing contribution to E 2 xp {t) as r — ► 0. Then we get 

a 2 ( T ) = e ~ rn ° V e^ 1 " (l - T ^\ + 0(r 2 ) = 

= cri(r) + rQ_ie- Tf20 V i e - rf7in + 0(r 2 ) , (13) 

where 0"i(r) is already obtained in (7). Using the known relation 

^i e - Qn = -ln(l-e- Q ) 

71=1 n 



-rn Q / t ,, rQ i 



we find 

a 2 (T)=a 1 (T)+T^ l v- r ^{hrr{h - -^ ) +0(7--) (H) 

So the regularized energy reads 



■cr p. . ' . ^ 1 . ^ 1 ^1 . ^1 . ^0 , ^n 



The contribution of the terms of order 0(l/n 2 ) (E* m ) is finite and can be found in any 
particular case. The divergent parts 

E quad (r) = 1— and E l ° 9 (r) = fi_! In r/i 

T ill 

is to be compensated using a certain renormalization scheme, which will be discussed later. 
Let us compare the result obtained above with another method of regularization — the 
^-function approach. In this case the energy is regularized as: 

El eta {s) = _^ Un -+- U ^s £ u -s (16) 

s— ►— 1 



where the arbitrary mass /i is included in order to restore the correct dimension. We assume 
that this mass equals /i in (16). For u^' = f2 x n + fi one has 

El eta {s) = - lim fi s+1 J2( Q i n + Q o)~ s -uj q = - lim Zl {s) - u , (17) 

n=l 

where 



zi(s) = // +1 ^(ftin + ft 



0; 



1 



= ?(fi^( 1 -^ + ^^feJ +0(l/n)j . (18) 

It is clear that the contribution of terms 0(l/n 3 ) vanishes as s — > — 1, so the regularized sum 
Ef eta (s) is determined by the leading three terms in the expansion (18). Then one finds: 



, +1 us + ij + 2 ^ +2 . 



.s+1 n-sA/„\ 6ii >•/ , i\ , "0 



*(*) = // +1 ^ S C(^) - 7^TTC(s + 1) + -^(a + l)C(a + 2) . (19) 



The values of ((z) analytically continued to the total real axis are known |J: 

«o) = -i, C(-D = -f = -i. (20) 

For ((z) in the vicinity of 1 we have 

lim C(z) = ~^-r + C , (21) 

z-*l z — 1 

where C = 0.5772156649... is the Euler constant. Thus taking the limit s — ► — 1 one obtains 

This expression reproduces explicitly the finite part of (9). The absence of the divergent 
term is due to the analytical continuation for C( z )- Taking into account the next term of the 
expansion in powers of 1/n (10) we find 

E* eta {s) = 
= — lim /i s+1 y^ ( Qin + fi H J — ^o = — lim ^2(5) — ^0 , (23) 

s-t— 1 ^ V 71 / »-*-! 



g a ( 8 )=z 1 ( 8 )- -^ C(^ + 2) (24) 



where 

sfi_i/i s+1 

(we write down only the terms which yield a non- vanishing contribution). The second term 
in (24) can be found using (21) and the expansion x £ = 1 + elnx + 0(e 2 ), e — s + 1, e — > 0. 
Then 

= n-.£s(-i + «)(i-^)(| + o) = -n-ig-i + o-h^) . (25) 



Therefore the regularized energy is of the form 



, _ fli no ng / fii \ /i 



£f"W = ^ + f + ^-«* + n-, lln^ + l l-n.^- + £7) (26) 



ilii^ oiiig uieii 1 1 _y jj/ 

_ = ^h ; i n7 = a (27) 



what coincides with (15) (except of the quadratic singularity E^ 10, ) provided that 

1 

T 



If one needs to take into account the contribution of terms of order O ( ^ J , then their sum 

(—E-f m ) should be added to (15) and (26). Therefore it is shown that the exponential and 
^-function regularizations provide the equivalent results for regularized Casimir energy in 

(1 + 1)1?. 

The singular part is: 

E dlv (r) = -^— + n_ 1 lnr f t. (28) 

It can be removed, for example, by means of the redefinition of the "bare" constants from the 
initial Lagrangian. It is very important to note that the way of renormalization is prescribed 
by the dependence of the extracted singularities on the geometric parameters (in our case the 
only geometric parameter is the bag's size R). We collect all singularities with the similar 
dependence from R, and then find the "contact term" to be redefined 0. 

Let us describe this procedure in the simplest case, i. e., MIT bag model with massive 
fermions in (1 + 1)D f!\, [10| . The Lagrangian reads 

C MIT = 9(\x\ < R)(iij>dil> - m F ijJip - B) + 0(\x\ > R)(iipdip - M F tptp) , (29) 

where B is the so-called bag constant which characterizes an excess of the energy density 
inside a hadron compared to the energy of the nonperturbative vacuum. Taking the limit 
M F — > oo in the exterior region |x| > R one gets the "bag", which is just the segment of 
real axis [-R, R\. The boundary conditions 

(±17* + l)ip{±R) = (30) 

lead to the spectrum 

' / 7T 7T x - 



Wn = \\ — n-\ ) +m%. (31) 

V V2i? ARJ f v ; 

Assume that the fermion mass vfip is small and expand the energy in powers of it. Then 
(with the accuracy up to m F ) one has 

^ (%rtl + i( fl ,fk) + ° K) ' (32) 

^! = 2R1 ^° = ar- ^ or n — 1 ^e expansion (4) can be found, where 

_ m% mLR , , 

n - = m = -i- (33) 



and for n = we have 



2 
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The divergent part of the energy ( according to (28)) is 

E div (r\R) = ( — \- + ^lnrJ 2R=2B'(r)R . (35) 

\ T 2 7T 2lX I 

In accordance with the approach of || [7| , the renormalization is performed by means of the 
redefinition of the bare bag constant B in the Lagrangian (29): 

B = B - B\t) . (36) 

The remaining parts of energy (15) are finite and can be found explicitly. The contribution 
of terms O (n~ 2 ) in the expansion (4) is determined by 

Ef in (R) = mp V { 1 — 

1 ' 2 ^ \VLrn + Q fiin 

m% ^-^ 1 2m%R , , . 

= ^rY—, r = — — (l-ln2). 37 

2fii ^i n(2n + 1) vr v ; v ; 

Therefore, the renormalized energy of the fermionic sea for (1 + 1)D massive MIT bag model 
as the function of the bag's size R and renormalized bag constant B reads (up to the terms 
of order mfp): 

E MIT (R) = 2B R-^+ T ^(l-2ln2 + \n-^-)+O(m 4 F ). (38) 

Taking the limit rnp — > (massless MIT bag model), one gets from (38) the well-known 
result of [p |. Note, that this configuration is unstable and tends to R — > 0. It is possible to 
make it stable by adding one valence fermion into the lowest level n = 0. Then the energy 
reads 

E MIT {R) = 2B R + HZL + HM U _ 2 In 2 + In ^-) + O(m^) . (39) 

3 Two-phase Hybrid Bag Model 

Now we are ready to study the hybrid bag model, in which the fermion ("quark") field 
interacts with scalar ("meson") field. Compared to the models considered in ||, here is no 
phase of the massless fermions. The Lagrangian reads 

C = iipdi/j + \ d^ip - 9(\x\ < R) (| m F [ij } e i75 ^] - B) - 

- 6{\x\ > R){V(tp) + \ M F [^,e^] )• (40) 

The commutator in the fermionic forms provides the charge-conjugation symmetry. For the 
sake of simplicity, we consider the scalar field <p(x, t) in the mean-field approximation (MFA) 



12], i. e., assume it to be a c— number function of the space-time variables; suppose also 
that it is independent of the temporal coordinate: tp(x,t) = tp(x). Besides this, we will 
take into account only the leading terms of the perturbative expansion in the chiral coupling 
constant [fL2 ]. In other words, the fermionic mass mp is assumed to be small inside the 
segment \x\ < R while the mass Mp in the exterior region |x| > R is infinitely large. Thus 
the fermion field at \x\ > R vanishes and the equations of motion for \x\ < R are 

(id-mpe il5V )^(x,t) = 0, (41) 






^,7 5 e i75 ^) sea , (42) 



where in r. h. s. of (42) the v. e. v. of the axial current is taken accordingly to MFA. For 
\x\ > R, the scalar field is determined by the non-linear equation 

- <p»(x) = V» . (43) 

The solution of the equations (41, 42) have been firstly found by Sveshnikov and Silaev and 
studied in detail in 0, so we will discuss it here only briefly. Let us assume that the scalar 
field is an odd soliton function and for \x\ > R it has the form: 

<p(x) = 7T (l - Ae~ mx ) , x > (44) 

where m is the meson mass, and (p(x) for x < is determined by the oddness. The boundary 
conditions are 

(±z 7 x + e^ Mx) )ij(±R) = , 

<p(±R±0) = tp(±R TO) , 

ip'{±R ± 0) = (p'{±R =p 0) . (45) 

The fermionic spectrum possesses the symmetry v — ► — v, where v 2 = k 2 + m 2 F . The 
corresponding unitary transformations of the wave function are x ~^ ijiX (here x — e n5(p ' 2 ijj), 
therefore the v. e. v. of the axial current in the r. h. s. of (42) is equal to zero and the 
solution of (42) for \x\ < R appears to be the linear function tp(x) = 2Xx (for detailed 
discussion of the self-consistent solution for this system see 0). The eigenvalues uo n can be 
obtained from the equation 

L _ c 2ikR m F + J{V + k) \ L _ ^ikR l V ~ k \ m F-J(v + k) \ = 

\ m F + %{y — k) J \ \v + k J m F — i{y — k) ) 

= (l- e -^ m F -i(u + k) \ L _ e „ 2ikR (v_-k\ m F + i{y + k) \ 
\ VfiF — i(u — k) J \ \u + k J nip + i[y — k) J 

where v = uo — A. Provided that the signatures of u n and v n are the same for all n, the 
fermionic Casimir energy is determined by (2) with the replacement u n — > v n . For the sake 
of being definite we regularize this expression by means of the exponential cutoff (3): 



rv n 



E^(r) = Jim ( - £ „ n e—) = lim ± £ e" 

V n >0 U ' U n >0 

The equation (46) can be written as 

m F sm2Rk + kcos2Rk = . (47) 

Using the expansion in powers of nip for k 

k = k + m F ki + m 2 F k 2 + 0(m F ) (48) 

one obtains 

u = k + m F k x + m 2 F I ~k 2 + -— - ) + 0{m F ) . (49) 

V 2fc / 



One can solve (47) in any order, what yields (up to O (m 3 )), 

71 n 2 1 lQRm 2 F o, . n . 

"" = 2l" + 4fi + ^ (1 + "" J! '(STT) - nH2„+ir + ° (m ' ) • (M) 

Expanding v n in powers of 1/n one gets 

i/ n = Q in + n + — + O ( ^-) , (51) 

where 

Sit = £-, Qo = ^=, fi-i = — (1 + m F i?) , (52) 

what differs from the massive MIT bag model by the term mpj-n in Q_ 1 which is independent 
of the bag's size R. For n = we have 

_ _ lQRm% , , 

i/o = fi + 2fi_i =-^ . 53 

7T ,:i 

Now one can use the formulae for the regularized energy (15). The finite fermionic energy 
should be of the form 

E F = E Cas + 2BR + A , (54) 

where B and A are the bag and cosmological constants, and Ec as is the divergent fermionic 
Casimir energy Hence the renormalization requires the redefinition of B (36) and A: 

A = A - — In Tfi . (55) 

71 

The contribution E^ m of the non-singular terms Of^j reads 

F fin = o v l + l6Rm F V X = 

- 1 y x n(2n + 1) + vr3 ^ (2n + 1)3 



16Rm 



fl_i2(l - In 2) + ^A , (56) 



7T" 

where A = 1.051799.... Therefore the finite energy of the fermionic sea is (provided that 
Ao = 0): 



E(R) = 2B R-^ + 



+ —(l + m F R) 

IX 



ln^^ + 1 -21n2 
2{iR 



16Rm2F -A . (57) 



7T 3 



A total energy of this system contains also the contribution of the scalar field. In order 
to find it we use the boundary condition for Lp and it's derivative (45) what yields 

2A = -^- . (58) 

mR + 1 y J 

By virtue of the virial theorem in the external region |x| > R one obtain | y? /2 (^) = V(ip), 
so the scalar field energy reads 

rR ( r-R r R\ tt-2™ 

E ip (R) = ± / dxip' 2 (x) + [ +/ )dx V '\x) = ——-. (59) 

J-R \J-oo Joo mR+l 



Note that the representation of the scalar field in the form (44) is valid only at the distances 
much larger than the soliton size, i. e. of order m _1 . This restriction allows us to postulate 
(in the framework of our model) the following relation between the meson mass m and the 
bag's radius: 

mR = K^ , K-n > 1 . 

Thus the contribution of the scalar field to the total energy of the bag is 

and the total energy consists of the fermionic and the bosonic parts: 

E tot (R) = E F {R) + E V {R) . (61) 

In terms of the dimensionless variables the total energy £ = E/irtp reads (x = mpR): 

S(x) = 2B lX - -?-+ 
48x 

ir (1 + K v )x 



1, 

-(1 + x) 



In - + l-21n2 

2x 



It is easy to see that this energy has the unique minimum and the configuration is stable. 

4 Conclusion 

The present approach to renormalization of the infinite ground state energy of the models 
with confined fermion fields in (1 + 1)1? is based on the analytical regularization of the di- 
vergent sums |13| (e. g., the exponential cutoff, or ^-function regularization), the extraction 
of the singular terms in the form of poles at the point of the physical value of regularization 
parameter, and subsequent absorbing of these singularities by means of redefinition of bare 
constants from the initial Lagrangian. In this framework the singularities are isolated un- 
ambiguously using the scheme analogous to MS in QFT. The remaining parts are finite and 
demonstrate the non-trivial dependence on the geometrical parameters of the model. The 
generalization of the proposed framework to the higher dimensions appears to be a separate 



problem and will be the subject of the future investigations |L4 . 
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